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Pion pair production in e + e annihilation 

Zhun LrQ and Ivan SchmidlQ 
Departamento de Fisica, Universidad Tecnica Federico Santa Maria, Casilla 110-V, Valparaiso, Chile 

We present an analysis of the process e + e~ — > 7* or Z* — > nn"/, in the kinematical region where 
y/s, the cm. energy of the e + e~ pair, is large but much below the Z-po\e. The subprocess 
7* or Z* — > tttvj can be described by the convolution of the hard scattering coefficient 7* or Z* — > 597 
and the general distribution amplitude of two pions qq —* tttv. In the case of neutral pion, the pro- 
duction through 7* is the dominant process, which can therefore be used to access the generalized 
distribution amplitudes (GDAs) of the pion, especially their C-even parts. The 7Z interference term 
provides an alternative approach to extract the weak mixing angle sin 9w through measuring the 
helicity asymmetry in the process e + e~ — > 7r°7r 7. In the case of charged pion pair production, the 
Bremsstrahlung process dominates and its interference with e + e~ — > 7* — > 7r + 7r _ 7 can be applied 
to study the process 7* — > nn~/ at the amplitude level. 

PACS numbers: 13.66.Bc, 13.88.+e, 14.40.Aq 



I. INTRODUCTION 

Generalized distribution amplitudes (GDAs) Q, Q 
are important nonperturbative objects for understanding 
how quarks and gluons form hadrons in hard exclusive 
processes. They describe the soft transition qq — > hh 
or gg — > hh, in the kinematical regime where the in- 
variant mass of the hadron pair hh is small compared 
to the hard scale of the process. The typical processes 
to which the GDAs formalism can be applied are the 
two photon process 7*7 — ► ^03^0, Q and hadron pair 
clctroproduction j*p — > p'hh jy, Fjj, for which factoriza- 
tion theorems have been proved [1] . In Ref . we gave 
an analysis of the process e + e~ — > Z — > irir-f at the Z- 
pole energy, with small pion pair invariant mass. The 
subprocess Z — > ttttj can be factorized into the convo- 
lution of a hard coefficient and the generalized distribu- 
tion amplitude of two pions (2-pion GDA) . In the case of 
charged pion pair production, the process is sensitive to 
the C-odd parts of 2-pion GDAs, and provide an oppor- 
tunity to access them in e + e~ annihilation. In this work 
we extend the analysis in Ref. Q to a generalized process 
e + e~ — ► 7* or Z* — ► 7T7T7, in the kinematical region where 
y/s is large but much below the Z-pole. Like the case of 
the production through Z boson analyzed in Ref. a 
similar factorization holds for the case of the production 
through a virtual photon, which gives the dominant con- 
tribution for neutral pion pair production in the kine- 
matical region we discuss here. Therefore, the investi- 
gation of the process e + e~ — > 7* — > 7r°7r°7 can provide 
detailed information about 2-pion GDAs, especially their 
C-even part. We also discuss the ^Z interference term, 
which contains the weak mixing angle sin 9w and can be 
used as an alternative approach to extract its value. The 
Bremsstrahlung process e + e~ — > 7*7 — > 7r + 7r _ 7 dom- 
inates in the case of charged pion production, and its 
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interference with e + e~ — > 7* 
tool for studying the process 7* 
level. 



7r + 7r 7 is a valuable 
7T7T7 at the amplitude 



II. ANALYSIS OF THE PROCESS 7* 



7T7T7 



We first study the subprocess of pion pair production 

7*(g) -> tt(pi) + ir(p 2 ) + 7(9') (1) 

in the framework of GDAs. It is convenient to investigate 
this subprocess in the pion pair center of mass frame 
(shown in Fig. QJ,), where the kinematics is expressed 
as (choosing the z axis along the momentum of the real 
photon) 
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respectively. Here v and v' are two lightlike vectors which 
satisfy: v 2 = v' 2 = 0, v ■ v' = 1, W is the invariant 
mass of the pion pair and q 2 = Q 2 = s, P is the total 
momentum of the pion pair. The skewness parameter £ 
is the momentum fraction of plus momentum carried by 
7r(pi) with respect to the pion pair: 



pj _ 1 + (3 cos e 



0=\/l- 



4ml 

7T 

W 2 



(3) 



In the kinematical region W 2 < s, the process given 
in Eq. QJ can be analyzed in the framework of GDAs, 
which has been applied in Ref. to the case where the 
virtual photon is replaced by the Z boson at y/s = Mz- 
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(b) 

FIG. 1: The kinematics of pion pair production, (a) the cm. 
frame of the pion pair, (b) the cm. frame of e + e~ . 

From Fig. [2^. and the one with the two photon vertices 
interchanged, the hadronic tensor of the process can be 
calculated as PHH 

T!?-^ = 1 1 d*xe- 4 «-*<7r( P i)7r( Pa )|TJ 1 S n (x)^ n (0)|0) 

(4) 

which has the same form as that of the two photonpro- 
cess 0,0, and C, W 2 ) is the GDA defined as Q: 

*> g (z,(,W 2 ) = J ^ e ^ p+ ^ 

x(7r(p 1 )7r(p 2 )|^(a ; -)7 + V'(0)|0}, (5) 

where z = k + / P + is the momentum fraction of the quark 
with respect to the hadronic system. 

As discussed in Refs. 0,0, a useful constraint on the 
GDAs comes from charge conjugation invariance of the 
strong interactions: 

$™(z, C, W 2 ) = -$7(1 - z, 1 - C, W 2 ). (6) 

Following this constraint one can define the C-even and 
C-odd parts (in j(| they are also called isoscalar and 
isovector parts respectively) of 2-pion GDA as Q 

$±(z, c, w 2 ) = \[%*{z, c, w 2 ) ± $J*(z, 1 - c, w^ 2 )], 

(7) 

where the superscript 7T7t represents 7r°7r° or 7r + 7r . Then 
$™(z, C, W 2 ) = *+(z, C, W 2 ) + <f>~(z, C, W 2 ). (8) 



The properties of &f(z,(^,W 2 ) under the interchanges 
z — > 1 — z and ^ — > 1 — can be easily derived from 
© 0: 

4>+(z,C,^ 2 ) = -$+(l-z,C,^ 2 ) 

= $+(z,l-(,W 2 ), (9) 
$-(z,C,^ 2 ) - $-(l-z,C,M/ 2 ) 

= -$-(z,i-c,w/ 2 ). (10) 

The £ — > 1 — C exchange corresponds to the interchange 
of the two pions, thus in the case in which the final two 
pions are 7r°7r° we get 

$f ff ° (z, C, W 2 ) = w ° (z, 1 - C, VF 2 ). (11) 

Therefore there is only a C-even part for (z, £, W 2 ). 

In the following we use $^(z, C W 2 ) to represent the C- 

even/odd part of $^ +w {z, C, H^ 2 ) respectively. Isospin 
invariance implies Q 

$f 7r °(z,C,^ 2 ) = $+(^C,VF 2 ). (12) 

Thus according to only <f>+(z, £, I/F 2 ) contributes to 
the process 7* — > 7T7T7. Here we only consider the contri- 
butions coming from the it and d quark. Another conse- 
quence of isospin invariance is 

*£ = *J, = (13) 

Eq. J3J shows that the process e + e~ — > 7* — > 7T7T7 can 
be also applied to probe <I> + , just like the two photon 
process 7*7 — > 7T7t, which has been proposed 0] to access 
Interest on the process e + e~ — ► 7* — » 7T7T7 relies 
also on the fact that this process is related to its crossed 
channel, the deeply virtual Compton scattering process 
(DVCS) en — ► e7T7 (shown in Fig. [2b), and is referred to 
as timelike DVCS. In the factorization of the DVCS pro- 
cess, the nonperturbative objects are the generalized par- 
ton distributions (GPDs) [13, EH, 

which arc important 
for understanding the parton angular momenta inside the 
nucleon [13 . The simplest DVCS process is eir — > e7T7. 
However it is difficult to use the pion as the target, which 
makes the pion GPDs difficult to be experimentally ac- 
cessed. It has been shown that there is a crossing 
relation that connects GPDs and GDAs. Therefore to 
measure the process e + e~ — > 7* — > 7T7T7, the crossed 
channel of the DVCS process, can also give experimental 
indication about the GPDs of the pion. In the kine- 
matical region W 2 <C s <C M|, the suitable facilities to 
investigate this process are B-factories. Finally, we re- 
mark that there is a suggestion n4| of using the transition 
distribution amplitude (TDA) [TJj formalism to describe 
the process e + e~ — > mvy. 
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(a) (b) 



FIG. 2: (a) Diagram of the process e + e nn~/, where V can be a photon or Z boson, (b) Diagram of the DVCS 

process eir — > eir^y, the crossing related channel of the process depicted in (a). 



III. PHENOMENOLOGY ABOUT THE 

e+e" -> 7T7T7 PROCESS 



vertex, with A e — ±1 the hclicity of the electron. The 
amplitude of V* — > 71717 can be calculated from 



A. Calculation of the cross-section 

In this section we study the phenomenology associated 
with the process 

e+(l) + e-(i') -> V*(q) -> 7r(pi) + tt( P2 ) + 7(9'), (14) 

where the intermediate vector boson V can be a photon 
or a Z boson, based on the analysis given in Section 
ITT1 and in Ref . Q . We will consider the helicity of the 
electron in the analysis. It is useful to define two Lorentz 
invariant variables: 



lq-q' 



y 



q-q r 



(15) 



The differential cross-section for the process e + e 
V* — > 77171 is expressed as 



d(J e + e — — r y * — >^y7T7T 



0W(s - W 2 ) 

64(2tt) 5 s 2 
xdWdn*dfl', 



1-4 



e+ e — > V * — *-7T7T7 I 



(16) 



in which (|p* |, f2*) is the momentum of pion 1 in the cm. 
frame of the pion pair, and fl' is the angle of the photon 
in the rest frame of the Z boson (that is, the cm. frame 
of e + e~). 

The amplitude of the process can be calculated from 



A e + e — *V* — >77T7T 



-A 



V* — >77T7T 



S-M*r l ' J 



(17) 



where = e7 A1 (l + A e 7 5 ) for lepton-7 vertex, and = 
e 7 M ( c \/ — 15 c a)(^ + ^75)/ sin2#w for lepton-Z boson 



Ai;^™ = e$e'( j >Tf(z, C, W\ (18) 



where e and e' arc the polarization vectors of the interme- 
diate virtual boson (photon or Z) and the real photon, 
respectively. In our reference frame these vectors have 
the form: 



(19) 



Using the notation of Ref. [T(| the cross-section of the 
process can be written as: 



da 



dW 2 d cos 8dcj)dy 32tts 4 



Y, L tXXa b - (20) 



a b 



The indices a, b can be 7 for the photon or Z for the Z 
boson. The relative propagator factors Xab are given by 

X 77 = !> ( 21 ) 

*r* = ^- sin2( 2^ )s _M r ( 22 ) 
Xzz = (x 7 z) 2 . (23) 



which correspond to the contributions from 7, 7Z inter- 
ference and Z boson. The lepton tensors are: 



L% = -2sC ab h{y)9T ~ isD ab I 2 {y)e^. (24) 
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FIG. 3: The 7 (solid line), -yZ interference term (dashed line) 
and Z (dotted line) contributions to the helicity-independent 
differential cross-section of the process e + e~ — > V* — > 7r°7r°7 
at W = 0.7 GeV, 9 = 10°. 



The notation used in the above equation is 



C 7Z 
Czz 

D Z z 
h(y) 
h{y) 



1, 

Cz 7 = (cy — c^A e ), 

(Cy + c/4) — 2CyC l A \ e , 



= x. 



D Z 7 = (CyA e - C l A ), 

(\-y + y 2 ), 

(l-2i/). 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 

(31) 
(32) 



The hadron tensors appearing in Eq. I|20|) are: 



-91 



J2e 2 q V q (cos0, W 2 ) 



(33) 



Kz = ~9T Re 



x 



9 

-ie , Im 



+ 



J2e q c q v V* (cos 6, W 2 ) 

Y,4V q (cos8,W 2 ) 

^2e q a v V*(co S 0,W 2 ) 
1 

^e^cosf?,^ 2 ) 

^e^V^cos^W 2 ) 

Y,e q c q A A q {cos9,W 2 ) 

Y,e q c v V; [cosO, W 2 ) 
^ e q c A A q (cos 0, W 



(34) 



-91 



L 9 



(35) 



where 



V g (cos6»,VF 2 



A 9 (cos6»,VF 2 ) 



, 2z-l 

z(l - z) 

dz ~, r 

z{l-z) 



<f>+(z,C,^ 2 ), (36) 



%{zX,W 2 ). (37) 



The vector and axial-vector couplings to the Z boson 
appearing in the above equations are given by: 



c v = T3 - 2Q l sin 2 9 W , 



(38) 
(39) 



where Q % denotes the charge and T\ is the weak isospin 
of particle i (i.e., T3 = +1/2 for i = u, and —1/2 for 
i = e~, d, s). 



B. Results for the neutral case 

In the previous subsection we have given the expres- 
sion for the cross-section of the general process e + e~ — > 
7*orZ* — > 7T7T7, in terms of 2-pion GDAs. There are 
interesting phenomenological implications coming from 
Eq. (|16|) . As demonstrated in Ref. 0, measuring the 
charged pion pair production in the kincmatical regime 
•y/s = Afz, where the main contribution comes from Z 
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FIG. 4: The W-dependent differential cross-section of the FIG. 5: The y/s dependent ratio of the helicity dependent 



process e T e 



7r°7r°7 at y/s = 10 GeV and 9 = 20° 



and helicity-independent cross-section of the process e e 
V* -> 7r°7r°7. 



boson, can provide an opportunity to access the C-odd 
part of 2-pion GDAs $~. In this work, however, we con- 
sider a different kinematical region W 2 <C s >C M§. 
Consequently a different phenomenology emerges: the 
production through 7* is dominant over the interfer- 
ence and Z boson terms, which is a result coming from 
the factors given in Eqs. J5TJ, and 

To avoid the competing contribution from the 
Brcmsstrahlung process, in this subsection we limit our- 
selves to the case of neutral pion pair production. There- 
fore the terms containing A q in the hadronic tensor l)24|) 
do not contribute to the cross-section, and the calcula- 
tion is simplified. Based on Q16|). we calculate the 7, 7Z 
interference and Z contributions to the unpolarized dif- 
ferential cross-section of the process e + e _ — ► 7* orZ* — > 
7r°7r°7 as functions of y/s. The results are shown in Fig. El 
by a solid line, dashed line and dotted line, respectively. 
For <!>+ , which is needed in the calculation, we adopt the 
simple model given in Q . The curves show that the con- 
tribution from the production through 7* is about several 
orders of magnitude larger than those through 7Z inter- 
ference and Z boson, in the regime W 2 <C s <C M§. 
Therefore it would be suitable to investigate <J>+ through 



ence term: 



the process e" 1 



7r 7r 7, in the regime where 



yfs is several GeV, where the cross-section is relatively 
large. 

The 7Z interference term is also interesting because it 
contains the weak mixing angle smOw- To isolate this 
term one needs to go beyond the helicity-independent 
cross-section, that is, to consider the helicity dependent 
cross-section of the process. The helicity dependent dif- 
ferential cross-section is dominated by the 7Z interfer- 



da(+X e ) — da(— A e ) 
2dW 2 dcos8 
a 3 /3(s - W 2 ) , 



6s 3 



-Xizc A X e 



x(e u c v + e d c v )(e 2 u + ej)\V u (cos9, W 2 )\ 2 , (40) 

since the Z boson term can be ignored compared to that 
of the "fZ interference term. 

As discussed previously, the unpolarized differential 
cross-section is dominated by the virtual photon produc- 
tion, which reads: 

da(+X e ) + da{-X e ) 



2dW 2 d(cos( 



a 3 P(s - 14^ 2 ) 
V2? 



(e 2 u + e 2 ) 2 \V u (co S 0,W 2 )\ 2 . (41) 



The VF-dependence of this cross-ection is shown in Fig. 01 
for y/s = 10 GeV and 9 = 20°. 

The ratio of the helicity dependent and helicity- 
independent cross-section, that is, the helicity asymme- 
try of the process is (A e = +1) 



R 



a(+X e ) -a(-\ e ) 
o-(+A e ) + ff(-Ae) 

_ l e uCy + P-dCy 
-2c A X~fZ 2 , n2 



10 



2 sin^ 9 



w 



&m 2 {29 w ) s-M z z 



(42) 



which depends on s and s'm 9w- Therefore the mea- 
surement on this helicity asymmetry provides an alter- 
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da/(dW 2 d(cose)d(|>dy)[fb GeV" 2 ] 




FIG. 6: Upper panel: the differential cross-sections of the Bremsstrahlung process, at i/s = 5 (solid line) and 10 GeV (dashed 
line). Lower panel: contributions to the differential cross-sections coming from the unpolarized interference term between the 
Bremsstrahlung process and the process e + e~ — * 7* — > 7r + 7r~7 at yfs — 5 (solid line), 10 GeV (dashed line) and helicity 
dependent interference term (dotted line) at yfs = 5 GeV. The first column: the W dependent differential cross-sections at 
9 — 20°, 4> = 180°, y — 0.3; the second column: the 9 dependent differential cross-sections at W — 0.7 GeV, <j> = 180°, y = 0.3; 
the third column: the (f> dependent differential cross-sections at W = 0.7 GeV, 9 = 20°, y = 0.3; The fourth column: the y 
dependent differential cross-sections at W = 0.7 GeV, 9 = 20°, <f> = 180°. 



native approach to extract the value of weak mixing an- cess has the form 
gle shi9w- In Fig. [S] we show the y/s dependent ratio 
R which is given in For y/s = 10 GeV, which is 

the typical cm. energy of CLEO or £?-factories, such ^ b a 3 3 (s — W 2 ) 

as BaBar and Belle, this ratio is about 0.8%, which is a = — ' \F 7T (W 2 )\ 2 



sizeable value. ' ' ' dW 2 d{cos6)d^dy 16ns 2 W 2 

a 

-4a; (1 - x) cos 2 



X((l-2s(l-a0) h ^ . sin 2 1 
2/(1 ^ V) 



Hv) 

V'! 

x(l -x)2 sin 2 6cos2(f>), (43) 



+ y/x(l-x)(l - 2x) 2[V) sin 26 cos 

vw - y) 



C. Results for the charged case 



In the case of charged pion pair production, besides the 
process e + e~ — > 7* — > 7r + 7r _ 7, the Bremsstrahlung pro- 
cess e + e~ — > 7*7 — » 7r + 7r~7 also needs to be considered. 
The cross-section of the former process is the same as that 
of neutral pion production, since only <1>J contributes. 
The differential cross-section of the Bremsstrahlung pro- 



where F 7T (W 2 ) is the timelike pion form factor. 

It is interesting to investigate the interference 0, Il7j 
between e + e~ — > 7* — > 7r + 7r~7 and the Bremsstrahlung 
process, because it can be applied to study the subpro- 
cess 7* — > 7r + 7r _ 7 at the amplitude level, similar to the 
case e7 — > e7T7r that has been discussed in Refs. 0,0|. 
Knowledge of the process 7* — > 7r + 7r _ 7 at the ampli- 
tude level can help to determine not only the magnitude 
of but also its phase. The form of the interference 
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term reads 



da 1 



dW 2 d{cos9)d(j)dy 



a 3 (3 2 {s-W 2 ) 



x( v / x(l-x)Re{F;v} cos6> 



- v) 

xRe{F*V} sin 9 cos 

h{y) 



-(l-x)Vy(i-y)] 



2^2/(1-2/) 



Im{F;y} sin0 sin (/.), (44) 



where V = ^2 q V q for q = u,d. This interference term 
can be obtained from taking the difference of the cross- 
section in the exchange (9, <j>) — ► (it — 9, tt + <fr). The last 
term in Eq. (|44|l is the hclicity dependent term, which can 
probe Im{F*V*}, and which together with the helicity- 
independent terms can be used to reconstruct the full 
complex amplitude of 7* — > 7r + 7r~7, in the region where 
Fir(W 2 ) is known. Conversely, these measurements could 
also be used in order to obtain information about the pion 
timclike form factor F 7r (W 2 ). 

In Fig. we give the curves of the differential cross- 
sections, coming from the Bremsstrahlung process and 
its interference with the process e + e~ — ► 7* — > it + 7t~j, 
respectively. For F^jW 2 ) we use the parametrization 
N = 1 given in Ref. We show the cross-sections for 
y/s = 5 (solid line) and 10 GeV (dashed line) in Fig. 
In the later case the terms containing 1 — x have sizeable 
contributions. In the lower panel of Fig. [5] we also show 
the contribution from the helicity dependent interference 
term (dotted line) for -/s = 5 GeV. 

Fig. shows that the contribution coming from the 
Bremsstrahlung process is 1 ~ 2 orders of magnitude 
larger than that from the interference term. Accord- 
ing to the shape of F„, the ^-dependent curves of 
both the Bremsstrahlung and interference contributions 
have peaks in the region W ~ m p . In the case of 9- 
dependence the curves of Bremsstrahlung and (unpo- 
larized) interference contributions increase with 9 from 
0° to 90°, and then decrease. The ^-dependent curve 
of the Bremsstrahlung contribution has a large value at 
6 = 180° and a small value at <h = 0°, while that of the 



(unpolarizcd) interference contribution has a large (ab- 
solute) value at (f> = 0° and a small (absolute) value at 
4> = 180°. The y-dependence of the interference contri- 
bution shows that the curves have maximum magnitude 
in the region y — ► or 1 (due to the factor \Jy{l — y) in 
the denominators of some terms in (|44[) h which is there- 
fore the suitable region in order to study the interfer- 
ence term. A similar situation occurs in the case of the 
Bremsstrahlung process. 



IV. SUMMARY 

We have given a detailed analysis of the process 
e + e~ — > 7* or Z* — ► 7T7T7, in the kincmatical region where 
\fs is large compared to the invariant mass of the pion 
pair but much below the Z-pole. A factorization simi- 
lar to which has been shown in Ref. [ij holds for the 7 
case: the amplitude of the process can be expressed as the 
convolution between the hard scattering 7* — > qqj and 
2-pion GDAs which describe the soft transition qq — > tttt. 
Thus the investigation on the process e + e _ — > 7* — > tttt"/ 
can provide detailed information about 2-pion GDAs, es- 
pecially their C-even part. The experimental measure- 
ment of this process can also provide experimental indica- 
tion of pion GPDs, by virtue of the fact that this process 
is the crossed channel of the DVCS process en — > eir-f, 
and there is a crossing relation between GPD and GDA. 
We studied the production of both neutral and charged 
pion. In the former case the 7Z interference term can 
be applied to extract the value of the weak mixing angle 
sin 9\y through the measurement of the helicity asymme- 
try in the process e + e~ — ► 7r°7r°7. In the last case the 
interference between the Bremsstrahlung process and the 
process e + e~ — > 7* — > 7r + 7r _ 7 is useful in order to study 
the subprocess 7* — > 7r + 7r~7 at the amplitude level. 
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